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Supercavitating  Flows:  Small  Perturbation  Theory 
and  Matched  Asymptotics 

Professor  K.V.  Rozhdestvensky 

St. Petersburg  State  Marine  Technical  University 
Lotsmanskaya  str.,  3 
190008.  St. Petersburg 
Russia 


Summary:  The  lecture  discusses  some  applications  of  the  theory  of  small  perturbations  as  applied  to  super- 
cavitating flows.  In  this  context  the  linear  theory  is  viewed  as  an  outer  expansion  of  a more  complete  nonlinear 
solution  of  the  flow  problem.  In  addition  to  comparing  linear  solutions  for  a supercavitating  flat  plate  for  differ- 
ent analogues  of  the  cavity  closure  models,  two  examples  are  considered  showing  how  to  account  for  the  presence 
of  local  flow  regions  where  the  perturbations  are  not  small.  In  the  first  example  a local  asymptotic  solution 
of  the  nonlinear  flow  problem  in  the  vicinity  of  the  leading  edge  is  matched  to  the  classical  linear  solution  to 
provide  a uniformly  valid  pressure  distribution  along  a supercavitating  flat  plate.  In  the  second  example,  the 
local  nonlinear  perturbation  of  the  otherwise  slightly  perturbed  flow  is  due  to  a spoiler  fitted  at  the  trailing 
edge  of  a flat  plate. 


1.  Introduction 


In  a nonlinear  steady  2-D  problem  formulation  for  a potential  flow  past  a supercavitating  body,  one  has  to  solve 
Laplace  equation  in  the  domain  occupied  by  the  fluid  with  the  following  conditions  on  the  flow  boundaries 


• Slip  condition  on  the  wetted  part  of  the  body 


&P  = &P\ 

dx  diA  &/’ 


V = Vb{x), 


(1) 


where  ip  is  perturbation  velocity  potential,  y = yi,(x)  is  a function,  describing  the  contour  of  the  wetted 
part  of  the  body  in  Cartesian  coordinate  system,  a;- axis  being  directed  downstream.  Note,  that  here  all 
quantities  and  functions  are  non-dimensionalized  with  use  of  the  characteristic  length  of  the  body  and  the 
velocity  of  the  oncoming  flow. 


• Dynamic  condition  on  the  cavity 

The  pressure  on  the  boundary  of  the  cavity  is  assumed  constant  wherefrom  the  corresponding  pressure 
coefficient  Cp  should  be  taken  equal  in  magnitude  and  opposite  in  sign  to  the  cavitation  number  a 

Cp  = CPc  = 1 - (1  + |^)2  - (J^y  = -a,  y = yc(x)  (2) 

where  y = yc(x)  describes  the  cavity  contour,  determined  in  the  course  of  the  problem  solution. 


• Kinematic  condition  on  the  cavity 

The  cavity  contour  is  a streamline,  and,  therefore,  it  should  be  subject  to  a slip  condition,  identical  to  (1). 

• At  the  infinity  the  perturbation  velocities  should  vanish. 

In  what  follows  we  first  consider  a steady  linearized  flow  problem  for  a supercavitating  foil1  using  different  linear 
analogues  of  the  cavity  closure  schemes.  Used  in  particular  are  the  analogues  of  the  closed  cavity  termination 
models  (Riaboushinsky  model,  Efros-Gilbarg  model,  Tulin  single-spiral  vortex  model),  as  well  as  those  of  Wu- 
Fabula  and  Tulin  double-spiral  vortex  termination  models.  Then  two  examples  are  presented  showing  how  the 
linearized  (outer)  solution  can  be  supplemented  by  a nonlinear  local  (inner)  solution  in  the  flow  regions  when  the 
perturbations  are  not  necessarily  small.  These  examples  include:  a uniformly  valid  solution  of  a flow  problem 
for  a flat  plate  at  an  angle  of  attack  (zero  cavitation  number)  and  of  that  for  a plate  with  a spoiler  at  the  trailing 
edge.  Both  examples  employ  the  method  of  matched  asymptotic  expansions  (MAE),  [2].  This  technique  consists 
in  finding  a local  (inner)  solution  in  the  appropriately  stretched  coordinates,  blending  it  smoothly  to  the  outer 
(linear)  solution,  and,  eventually  obtaining  a uniformly  valid  solution  by  additive  composition  of  the  inner  and 
outer  solutions. 


1The  term  “supercavitating” , as  understood  here,  implies  that  the  cavity  extends  beyond  the  trailing  edge  of  the  foil 
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2.  Linear  solution  of  the  problem  of  a supercavitating  flow  past  a thin  foil  with 
different  closure  schemes 


Consider  a linear  problem  of  a steady  supercavitating  flow  past  a flat  plate  with  different  closure  schemes.  In 
what  follows  all  quantities  and  functions  are  rendered  nondimensional  with  use  of  the  chord  c of  the  plate  and 
the  velocity  of  the  incoming  flow  Ua.  To  avoid  complicated  derivations  when  explaining  the  essential  points  of 
the  lecture,  one  assumes  that  the  cavity  detaches  from  a sharp  leading  edge  of  the  supercavitating  foil  of  zero 
thickness.  In  most  of  the  examples,  discussed  herein,  the  foil  is  represented  by  a flat  plate  at  at  an  angle  of 
attack. 

Assume  that  the  flow  perturbations  are  small.  Expanding  previous  nonlinear  formulation,  one  easily  shows 
that  the  linear  flow  problem  for  the  perturbation  velocity  potential  ip  = <p(x,  y)  is  governed  by  the  following 
equations 


• Laplace  equation 


d2ip  d2ip 
d^  + dy2 


( x , y)  € T>; 


• Flow  tangency  condition  on  the  wetted  part  of  the  foil 


dip  = d yQ 
dy  d* 


y = 0 — 0,  x € (0, 1), 


or,  for  a flat  plate  y0(x)  = —ax, 


dip 


y = 0 — 0,  are  (0,1), 


where  a is  angle  of  attack  (in  radians); 

• On  the  boundary  of  the  cavity 

As  the  linearized  pressure  coefficient  is  approximately  equal  to 


(3) 

(4) 


r,  _ 2 (p-Po)  _ 0d<p 

PU2  ~ zdx ■ 


and  the  cavitation  number  a is  defined  as 


_ 2 {pc  po ) 

pu2„ 


(5) 

(6) 


where  pc  and  p0  are  pressures  in  the  cavity  and  at  the  upstream  infinity.  Then,  the  perturbed  horizontal 
velocity  on  the  boundary  of  the  cavity  should  be 


dtp  1 
dx  ~ 2 <T 


y = 0 ± 0,  xe(0,l)  and  x6(l,l), 


(7) 


where  / is  the  cavity  length  measured  from  the  leading  edge  of  the  foil  and  related  to  the  chord  of  the  foil; 

• Condition  at  infinity 

V<^  = *§r  ->0,  f°r  x‘2  + y2^°°i  (8) 

where  i and  j are  unit  vectors  of  the  axes  x and  y correspondingly. 

It  should  be  noted  that  the  linearization  implies  that  both  a and  a tend  to  zero  at  the  same  speed,  i.e.  a = O(a). 
In  order  to  use  powerful  methods  of  the  complex  functions’  theory  it  is  convenient  to  re-write  the  flow  problem 
formulation,  described  previously,  in  terms  of  complex  variables  and  functions.  Instead  of  the  perturbation 
velocity  potential  ip  = ip(x,  y),  introduce  a complex  potential  F = F{z)  = ip(x.y)  + i tp(x,y),  where  i = y/—l 
is  an  imaginary  unit,  z = x + i y and  tp{x,  y)  is  a stream  function.  In  this  case  the  solution  of  the  problem  is 
reduced  to  finding  an  analytic  function  AF/Az  = w(z)  = u(x , y)  — iv(x,  y)  in  the  domain  V of  the  flow  with  the 
following  boundary  conditions  on  the  flow  boundaries 

• On  the  wetted  part  of  the  plate 


^w(z)  = —v(x,  y)  = a,  y = 0 — 0,  x 6 (0, 1); 


(9) 
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On  the  cavity 


• At  infinity 


5ftw(z)  = u(x,  y)  = -a,  y = 0 ± 0,  x 6 (0,  l ) and  x € (1, I); 


w(z)  — > 0,  for  z — > oo. 


£ 1.  Cavity  closure  models  with  a square  root  singularity 

In  the  nonlinear  formulation  this  type  of  the  cavity  closure  models  is  characterized  by  presence  of  a stagnation 
point  in  the  region  of  cavity  termination,  i.e.  that  of  Riaboushinsky  and  Efros-Gilbarg  (see  Fig.  1(1)  and  Fig. 
1 (2)).  As  discussed  in  [1]  a linear  analogue  for  all  such  models  contains  a square  root  singularity  near  the  point 
of  the  cavity  termination.  Consider  a corresponding  flow  problem  solution  following  [1]. 

The  linearized  “physical”  complex  plane  z of  the  flow  past  a foil  with  a trailing  cavity  of  finite  length  is  illustrated 
in  Fig.  2.  Note  that  in  this  plane  the  plate  plus  the  cavity  are  represented  by  a slit  y = 0 ± 0,  x 6 [0,  l) . To  solve 
the  problem,  map  the  exterior  of  the  slit  in  the  complex  plane  z onto  the  upper  half  of  the  auxiliary  complex 
plane  £ = £ + ir/  (see  Fig.  2)  with  help  of  the  function 


C = 


a = \/T— T. 


The  inverse  function,  i.e.  z = z(£),  can  be  written  as 

K2 

The  correspondence  of  the  points  in  z and  £ planes  can  be  seen  from  Fig.  2.  Note  that  the  point  z = l.  which 
represents  the  closure  point  of  the  cavity  in  the  “physical”  plane,  passes  into  a point  £ = oo  in  the  auxiliary 
halfplane  3£  = £ < 0.  In  fact,  for  z — ► l 

C « 00 • (14) 


On  the  other  hand,  the  infinity  of  the  “physical”  plane,  i.e.  z = oo  passes  into  the  point  £ = — ia.  Eventually, 
the  wetted  part  of  the  slip  coincides  with  the  interval  £ € (0, 1)  whereas  the  upper  and  the  lower  boundary  of 
the  cavity  in  z plane  are  mapped  onto  the  negative  rj  = 0 — 0,  £ < 0 and  positive  y = 0 — 0,  £ 6 (1,  oo)  parts  of 
the  real  axis  in  the  auxiliary  plane. 

In  the  context  of  the  closure  models  it  is  important  to  specify  the  anticipated  behavior  of  the  complex  conjugate 
velocity  w(Q  in  the  vicinity  of  the  cavity  closure  point.  It  can  be  shown  that  the  linear  analogues  of  the  nonlinear 
closure  models,  containing  a stagnation  point  within  the  cavity  termination  zone,  are  characterized  by  a square 
root  singularity  of  the  perturbed  velocity  at  the  corresponding  closure  point  (i.e.  at  z = 1).  Assuming  that 

for  z — » l w(z)  « — t==.  (15) 


Accounting  for  the  expression  (14),  it  means  that  in  the  atixiliary  plane 

w[z(0]^-^  = iAC,  (1G) 

i avl 

where  A is  a real  constant  to  be  determined. 

It  is  convenient  to  solve  the  boundary  problem  for  the  function 

w*{z)  = w(z)  - (17) 

For  this  function  5ft w*(z)  = 0 on  the  boundary  of  the  cavity,  i.e.  on  y = 0 — 0,  £ < 0 and  1 < £ < oo.  The 
latter  fact  simplifies  the  final  expression  for  the  solution.  Note  that  the  boundary  conditions  for  w*(Q  on  the 
real  axis  of  the  auxiliary  complex  plane  £ are  mixed,  prescribing  3w*(£)  on  one  part  of  the  axis  and  5 Rw*(Q  on 
another  part  of  the  axis,  see  Fig.  2. 

Assuming  a quarter-root  singularity  w(z)  = Q{z~1^)  at  the  leading  edge  and  square  root  zero  w(z)  = 0[(1  — 
z)1/2]  at  the  trailing  edge,  and  using  Keldysh-Sedov  formula,  one  can  derive  the  solution  for  the  complex 
conjugate  velocity  in  the  form 


, . 1 a / 1 — C 


i-£C-£ 


d£  + ih 


+A^ca-o- 
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Note  that  B and  A are  real  constants  to  be  determined,  and  the  last  term  in  the  formula  is  a permissible  solution 
accounting  for  the  required  behaviour  of  the  function  w(Q  for  £ — > L4£. 

After  integration,  one  can  obtain 

w(0=\cr  + (B-a)J^  + A^/C{r^0.  (19) 

The  solution  thus  obtained  contains  three  unknown  (real)  parameters  a,  B and  A.  The  first  relationship  between 
these  parameters  can  be  obtained  by  applying  the  infinity  condition  (8).  For  £ — > — ia  the  conjugate  complex 
velocity  should  equal  zero,  i.e. 

w(—ia)  = 0,  (20) 

The  latter  complex  equation  is  equivalent  to  the  two  real  equations 


5ftw(— ia)  = 0,  9w(— ia)  = 0. 


(21) 


The  remaining  relationship  necessary  to  determine  all  three  unknowns  of  the  solution  can  be  obtained  from 
the  closure  condition,  which,  in  other  words,  is  a requirement  that  the  cavity  should  have  finite  length.  The 
corresponding  equation  can  be  written  as 


/ Sr[w(z)]  d z = Q j)  w{z)  dz 


$^w(£)^d£  = 0, 

£C 


(22) 


where  contour  Cz  encloses  the  slit  j/  = 0±0,ie  (0,  /)  and  is  passed  in  the  clockwise  direction,  and  the  contour 
encloses  the  point  £ = — ia  and  is  passed  in  the  clockwise  direction.  The  derivative  dz/d£  can  be  derived 
from  the  expression  (13)  in  the  form 


dz  _ 2/£  2l(3 

d£  £2  + ®2  (C2  + °'2)2 


(23) 


Treating  the  equation  (22)  is  reduced  to  calculation  of  the  imaginary  part  of  the  residue  of  the  function 
w(£)dz/d£  at  the  point  £ = — ia,  i.e. 


9^'W(C)|dC  = a{27riRes[wC)|]|c=_ja}.  (24) 


Calculating  the  residue  with  use  of  the  Taylor  series  of  w(( ) in  the  vicinity  of  £ = — ia 

= (C  + ia)+0[(C  + ia)2]  (25) 

and  the  expression  for  the  derivative  of  the  inverse  mapping  function  2 = z(£),  one  comes  to  the  following  form 
of  the  equation  (22) 

MfLJ-0- 

The  lift  coefficient  can  be  derived  by  integrating  the  pressure  over  the  slit  in  the  “physical”  plane,  so  that 

cy  = = m / w(-z)  iz  = m / w(c)  (if  dc-  (27) 

cz  c( 


Similarly  to  the  previous  derivations  leading  to  the  equation  (22),  the  calculation  of  the  lift  coefficient  is  reduced 
to  that  of  the  real  part  of  the  residue  of  the  function  w(£)dz/d£,  which  yields 

(28) 

The  results  obtained  for  the  case  of  a flat  plate  at  an  angle  of  attack  a have  quite  a simple  form.  After  some 
algebra  the  expression  for  the  lift  coefficient  becomes 

Cy  = *al{-^L-l),  (29) 

The  relationship  between  cavitation  number  and  length  of  the  cavity  is  given  by  the  formula 

_ a _ 2 _ 2 

a a y/l  — 1 


(30) 
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Note  that  both  of  the  proceeding  formulae  imply  that  l > 1.  Using  the  expressions,  describing  parametric 
dependences  of  the  lift  coefficient  and  the  ratio  a /a  on  the  nondiinensional  length  of  the  cavity  l,  one  can 
obtain  the  following  formula 


c„ 


4+1 

a2 


2\/<t  2 \ 


(31) 


2.2.  Wu-Fabula  closure  model 


The  solution,  corresponding  to  Wu-Fabula  cavity  closure  model  (see  Fig.  1(4)),  can  be  easily  derived  from  the 
previous  solution  by  requiring  a “smooth”  termination  at  the  point  z = l,  or  ( = oo.  This  is  achieved  by  putting 
the  real  constant  A equal  to  zero.  The  corresponding  linearized  “physical”  plane  is  shown  in  Fig.  3.  The 
resulting  expressions  for  the  lift  coefficient  and  the  ratio  d = <j/a  versus  nondimensional  length  of  the  cavity 
are  given  by  the  formulae 


Cy_  _ ir\Tl 

a y/l+y/l-  1 ’ 
_ (T  _ 2 

o y/J  -f-  y/l  — 1 


(32) 

(33) 


2.3.  Tulin  double  spiral-vortex  cavity  closure  model 

In  the  double  spiral-vortex  model  introduced  by  Tulin  (Fig.  1(3))  it  is  assumed  that  the  cavity  termination  is 
followed  by  a wake  extending  to  the  downstream  infinity,  [1],  Therewith,  the  cavity  pressure  turns  abruptly 
into  that  of  unperturbed  flow.  The  corresponding  linearized  flow  problem  formulation  is  shown  in  Fig.  4.  Note, 
that  the  linearization  implies  that  the  abscissas  of  the  cavity  termination  points  on  the  upper  and  the  lower 
“banks”  of  the  cut  are  identical.  To  formulate  the  relevant  boundary  problem  for  the  perturbed  conjugate 
complex  velocity,  map  the  exterior  of  the  semi-infinite  cut  in  the  “physical”  plane  z onto  the  lower  half  plane 
= £ < 0 of  the  auxiliary  plane  £.  This  can  be  done  by  means  of  the  function 

C = ~\fz  (34) 


where  the  selected  branch  of  the  root  transfers  the  point  z = 1 + i(0  — 0)  into  the  point  £ = — 1 + i(0  + 0).  The 
correspondence  of  the  points  in  z and  £ complex  planes  is  indicated  in  Fig.  4.  Note,  that  the  termination  points 
of  the  cavity  z = l + i(0  =F  0)  are  mapped  respeetingly  into  the  points  £ = ±6,  where  b = \fl.  The  boundary 
conditions  for  the  conjugate  complex  velocity  'tc(C)  on  the  real  axis  = £ of  the  auxiliary  plane  are  shown 
in  Fig.  4.  Using  Keldysh-Sedov  formulae  to  solve  this  mixed  boundary  problem,  one  can  derive  the  following 
expression 


a 
2 n 


■«> -‘if-r 


-zc-z 


1-C 


<f*/V 


J d^_ 

Z-1C-Z 


+ B 


1-C 


(35) 


where  B is  a real  constant  to  be  determined.  Integrating,  one  obtains 


w(  0 


(B  — a)  * 


+ ia, 


(36) 


where 


T(C,Z)  = 


,/ZS 


_Z dC_ 

z-n-z 


= -2  111 


V2 


+ 


x/clc^ry  + v^ic^Ty 


(37) 


Requiring  that  the  perturbations  vanish  at  the  infinity,  one  obtains  B = 0.  After  some  additional  algebra  the 
expression  for  w(Q  acquires  the  following  form 


\Jb  + y/b  — 1 

Vb  + y/b+1 


ln^  + 2in^H)  + ^l 

C-b  v^7TI)  + y/6(C3T)J 


(38) 
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Finally,  the  expression  for  the  lift  coefficient  is  obtained  as 

— - = aby/b(y/b  + 1 — y/b  — 1),  a = — (39) 

a ' a 

where  b = yfl.  The  ratio  a versus  I is  described  by  the  following  formula 

- = - = - \Vb(Vb^T  + VbTT)  in ^+v^±Il  (40) 

O O 7T  L ' Vb+yfb^li 

The  previous  two  expressions  can  be  viewed  as  a parametric  dependence  of  the  lift  coefficient  upon  cavitation 
number  and  angle  of  attack,  or  rather  as  Cy/a  as  a function  of  a /a  for  the  linear  supercavitating  flow  past  a 
fiat  plate  with  a Tulin  double  spiral-vortex  cavity  closure  scheme.  Note  that  in  order  to  find  the  relationship 
between  the  cavitation  number  a and  the  length  l of  the  cavity,  it  was  assumed  herein  that  the  total  drag  of 
the  cavity+wake  system  should  be  equal  to  zero.  The  latter  statement  can  be  formally  expressed  as 

9 / w2(C)^dC  = °-  (41) 
<=oo 

The  corresponding  calculation  is  reduced  to  finding  the  resudue  of  the  integrand  at  the  infinity  in  the  auxiliary 
plane 

Plotted  in  Fig.  5 are  curves  of  lift  coefficient  related  to  angle  of  attack  versus  cavitation  number  related  to 
angle  of  attack,  i.e.  Cv/a  = f(a/a)  for  different  linear  analogues  of  the  cavity  closure  scheme  for  the  case  of  a 
flat  plate.  Presented  in  the  Figure  are  the  results: 

1.  For  Tulin  open  closure  scheme  featuring  two  double-spiral  vortices, 

2.  For  a cavity  scheme  with  a square-root  singularity  of  perturbation  velocity  at  the  cavity  termination 
points,  and 

3.  Wu-Fabula  cavity  closure  scheme  (no  singularity  at  the  cavity  termination  point) 

Plotted  in  the  same  Figure  are  experimental  points,  obtained  by  F.F.  Bolotin  and  E.B.  Anoufriev  for  a series 
of  segment  foils  with  different  relative  thickness  4.2,5  and  6%,  [3].  It  should  be  noted  that  in  the  case  of  the 
developed  cavitation  starting  from  the  leading  edge,  the  flow  past  the  segment  foil  is  equivalent  to  that  of  the 
flat  plate.  It  follows  from  the  comparison  of  test  data  with  the  calculated  results  that  the  correlation  is  fair  not 
only  for  the  case  of  developed  cavitation  (long  cavities),  but  also  for  the  transitional  regime  when  1 < l < 1.5. 
Tulin  double-spiral  vortex  scheme  provides  satisfactory  results  for  the  lift  coefficient  of  flat  plate  up  to  a /a  « 6, 
i.e.  practically  up  the  boundary  of  supercavitation.  The  experimental  points  are  seen  to  be  located  between 
the  schemes  1)  and  2).  Wu-Fabula  scheme  allows  to  obtain  the  magnitudes  of  Cy/a  only  for  a/a  < 2,  i.e.  for 
sufficiently  long  cavities.  For  a/a  close  to  2,  Wu-Fabula  scheme  gives  somewhat  excessive  magnitudes  of  the 
lift  coefficient. 

2.4 ■ The  case  of  zero  cavitation  number  - analogy  with  a fully  wetted  foil. 

The  simplest  albeit  practical  case  corresponds  to  zero  cavitation  number.  Therewith  the  cavity  becomes  semi- 
infinite. The  “physical”  flow  domain  in  z-plane  is  transformed  onto  an  auxiliary  lower  half-plane  = rj  < 0 
with  the  mapping  function  £ = — </z,  discussed  previously.  As  before,  one  can  easily  construct  the  solution  of 
the  ensuing  boundary  problem  for  w(Q  with  use  of  the  Keldysh-Sedov  formula.  Restricting  the  analysis  to  the 
case  of  a flat  plate,  one  comes  to  the  following  simple  formula 


which  satisfies  the  requirement  that  the  perturbed  velocity  should  vanish  at  the  infinity.  The  latter  property  of 
(42)  can  be  easily  verified.  In  fact  for  large  z — > oo,  C — ► oo 


In  particular,  the  perturbed  velocity  distribution  on  the  wetted  (lower)  side  of  the  plate  (£—►£  — i0,£  G (0, 1)) 
becomes 


m(|,  0 — 0)  = —a 


(44) 
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Referring  to  the  classical  theory  of  a thin  foil,  one  can  conclude  that  the  above  result  coincides  with  the  perturbed 
velocity  distribution  on  the  lower  side  of  a flat  plate  placed  in  a uniform  flow  in  (-plane.  This  suggests  a concept 
of  an  equivalence,  existing  between  the  original  supercavitating  flow  in  z plane  and  a fictitious  non-cavitating 
(fully  wetted)  flow  in  (-plane,  which  was  first  indicated  by  Tulin. 

It  should  be  noted  that  since  the  fluid  domain  has  been  mapped  entirely  into  the  lower  part  of  the  (-plane,  the 
region  9(  = 77  > 0 has  no  physical  significance.  Mathematically,  it  can  be  regarded  as  the  second  Riemann 
sheet  of  the  0-plane,  or  as  the  domain  of  the  cavity.  Owing  to  the  fact  that  on  the  boundary  of  the  cavity  in  ( 
plane  the  real  part  of  the  conjugate  complex  velocity  w(()  is  zero,  one  can  use  the  Schwartz  reflection  principle 
to  continue  w(()  across  the  real  axis  5ft(  = ( (excluding  the  cut  ( 6 (0, 1)),  provided  u is  an  odd  function  of  77. 
To  ensure  that  w(()  = u((,  r/)  + i'f((,  77)  is  analytic  in  the  entire  (-plane,  it  follows  from  the  Cauchy-Riemann 
equations  that  v should  be  an  even  function  of  77.  Thus,  the  appropriate  reflections  into  the  upper  half-plane  of 
( are  given  by 

ufarj)  = -u((,  -77),  (45) 

v(&  v)  = 7 ;(&  -77)  (46) 

The  reflected  boundary-value  problem  shown  is  mathematically  identical  to  the  lifting  problem  of  a fully  wetted 
thin  foil.  In  both  cases  a Kutta-Zhukovsky  condition  should  be  imposed  at  the  trailing  edge,  and  the  perturbation 
velocity  should  vanish  at  the  infinity.  Thus,  the  solution  of  the  supercavitating  flow  problem  (a  = 0)  in  the 
complex  plane  z 'written  in  auxiliary  complex  variable  ( is  identical  to  that  of  a fully  wetted  foil  in  ( ] -plane . This 
analogy  enables  to  obtain  the  coefficients  of  hydrodynamic  forces  and  moments  of  the  supercavitating  foil  from 
those  of  a fictitious  fully  wetted  foil.  For  example,  to  obtain  the  lift  coefficient  Cv  of  the  supecavitating  flat 
plate,  one  has  to  integrate  the  pressure  coefficient  along  the  wetted  part  of  the  plate,  i.e. 


°v  pU2c 


o V f1  f1 

2 = / Cp(x,  0 — 0)  dx  = — 2 / u(x,  0^0)  Ax. 
pU  c ,/0  ,/0 

2', 


Passing  over  to  the  integration  in  (-plane  (x  = (2).  one  derives  from  the  previous  line 

Cy  = - 2 [ u(x,  0 - 0)  dx  = -2  / «,(£,  0 - 0)  ~ d(  ==  -4  [ «,((,  0 - 0)£  d( 
Jo  Jo  Jo 


(47) 


(48) 


Accounting  for  the  adopted  reflection  of  u into  the  upper  half- plane  C,  one  can  re-write  (48)  in  the  following 
way 

cv  = -4  f «(!,  0 - 0)(  d(  = 2 / K(.  0 + 0)  - 7 /,((,  0 - 0)](  d(  = 

f [Cp((,  0 - 0)  - Cp(t  0 + 0)](d(  = , (49) 

Jo 

where  Cm(  is  a coefficient  of  the  longitudinal  hydrodynamic  moment  of  a fictitious  fully  wetted  foil  of  length 
Cf  = \fc  in  (-plane.  This  coefficient,  calculated  with  respect  to  the  leading  edge  of  the  fictitious  plate,  is  defined 

_ 2 Mc  _ 2M< 

Gm<  " PU2C2  - PlP  y/C\fc 

Similarly,  one  can  show  that  the  coefficient  of  the  longitudinal  hydrodynamic  moment  of  a supercavitating  foil 
(cr  = 0)  is  equal  to  the  coefficient  of  the  third  moment  of  the  fictitious  fully  wetted  foil. 


Cm=  f Cp(x,  0 — 0)x  Ax  = [Cp((,0-0)-Cp((,0  + 0)](3 


d( 


It  can  also  be  shown  that  the  drag  coefficient  for  the  supercavitating  foil  can  be  found  as 

Cx  = 


(51) 


(52) 


This  useful  equivalence  was  first  established  by  Tulin  and  Burkart  in  1955.  For  a supercavitating  flat  plate 
(cr  = 0)  at  a given  angle  of  attack  a the  corresponding  calculations  give 


_ 7 ra 
Cy~  2 ’ 


C„,  = 


5770 

"32“  ’ 


CT  = 


(53) 


It  is  interesting  that  the  lift  coefficient  for  the  supercavitating  flat  plate  is  four  times  less  than  that  of  the  fully 
wetted  plate  and  two  times  less  that  for  the  case  of  a flat  plate  gliding  on  the  surface.  In  comparison  to  the 
fully  wetted  flat  plate,  here  the  centre  of  pressure  is  shifted  from  the  quarter-chord  point  to  a position  5/16 
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of  the  chord  downstream  from  the  leading  edge.  Note  that  the  drag  coefficient  of  the  supercavitating  plate  is 
seen  to  be  equal  to  Cya,  i.e.  the  total  force  is  normal  to  the  plate.  This  is  different  from  the  fully  wetted  case, 
where,  according  to  D’Alembert  paradox,  the  vector  of  the  force,  acting  upon  the  foil,  is  strictly  vertical  due 
to  the  suction  force  at  the  leading  edge  and  the  resulting  drag  is  zero.  In  the  supercavitating  case  there  is  no 
suction  force,  since  the  square  root  singularity  of  the  equivalent  fully-wetted  flow  in  the  (-plane  is  reduced  to  a 
quarter-root  singularity  in  the  physical  z-plane. 


3.  Uniformly  valid  asymptotic  solution  of  the  problem  for  the  flow  past  a super- 
cavitating flat  plate 


The  linear  theory  reveals  a quarter-root  singularity  of  the  perturbation  velocity  at  the  leading  edge  of  a super- 
cavitating foil.  This  deficiency  can  be  corrected  by  a special  consideration  of  the  local  region  of  the  flow  near 
the  leading  edge,  as  first  done  in  [4]. 

3.1.  Linear  solution  as  a leading  order  outer  expansion 

Consider  the  perturbation  velocity  distribution  along  the  wetted  (lower)  side  of  the  supercavitating  (a  = 0)  flat 
plate,  following  from  (44).  In  the  auxiliary  variable  5R(  = £ G (0, 1)  the  horizontal  component  of  w equals  to 


u(£,  0 - 0)  = 5IMC) 


(54) 


As  follows  from  the  mapping  function  ( = —\fz  the  real  variables  x and  £ are  related  to  each  other  as 

£ = for  y = 0 ± 0 


(55) 


so  that  on  the  wetted  side  of  the  supercavitating  plate  x G (0, 1),  j/  = 0 — 0 the  expression  (54)  can  be  re-written 
in  terms  of  “physical”  x coordinate 


u(x,  0 — 0) 


(50) 


This  expression  shows  explicitly  the  quarter-root  singularity  of  the  perturbation  velocity  (pressure)  in  the 
vicinity  of  the  leading  edge  of  the  wetted  side  of  the  supercavitating  flat  plate.  In  fact,  for  x — ► 0 + 0 


u(x,  0—0) 


— x Cp(x,  0 — 0)  ~ 


-zv-r-1/ 4 


(57) 


Considering  this  linear  theory  solution  as  an  outer  expansion  of  a complete  solution  of  a nonlinear  problem  for 
a flow  past  a supercavitating  flat  plate  for  a — > 0,  i.e. 


u(x,  0 — 0)  = u °(x,a)  = au°(x ) + O (a2),  u°(x)  = 0(1) 


(58) 


one  can  note  that  this  asymptotic  expansion  loses  uniform  validity  for  sufficiently  small  x when  the  product 
au°(x ) acquires  the  order  of  0(1).  This  takes  place  at  distances  from  the  leading  edge  of  the  order  of  x = 0(a4). 
For  the  purpose  of  further  analysis,  one  would  need  to  have  the  expression  for  the  ordinates  of  the  upper 
boundary  of  the  cavity  from  the  linear  theory.  Recalling  that  for  2 = x + i(0  + 0),  x > 0 the  £ and  x coordinates 
are  related  as  £ = — y/x,  £ < 0,  one  arrives  at  the  following  outer  expression  for  the  vertical  component  of  the 
conjugate  complex  perturbation  velocity 


v(x,  0 + 0)  = c3w(x)  = v°(x,  a) 


—a  + avi(x,0  + Q)+O(a2),  v°  = 0(1)  (59) 

The  outer  expansion  for  the  ordinates  of  the  upper  boundary  of  the  cavity,  measured  with  respect  to  the  plate 
can  be  determined  with  use  of  the  following  formula 

y°(x,  0 + 0)  = a f v°(x,  0 + 0)  d.r  + 0(«2)  = ay°(x,  0 + 0)  + 0(a2).  (60) 

Jo 


Integrating  this  expression  with  account  of  the  formula  for  v°(x,  0 + 0)  gives  the  following  expression  for  the 
outer  contour  of  the  upper  boundary  of  the  cavity 


y°(x,0  + 0)  ~ m/J(a:,0  + 0) 


(61) 
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3.2.  Inner  expansion  in  the  vicinity  of  the  leading  edge 

As  follows  from  the  proceeding  estimates,  the  linear  solution  is  not  valid  in  the  vicinity  of  the  leading  edge  of  the 
supercavitating  plate,  having  dimensions  of  the  order  of  0(a4).  As  a result  of  this  nonuniformity  the  pressure 
on  the  wetted  side  at  the  leading  edge  becomes  unrealistically  infinite.  This  also  means  that  within  the  initial 
characteristic  scale  of  the  order  of  the  length  of  the  plate  the  stagnation  point,  located  close  to  the  leading  edge 
is  invisible.  In  order  to  the  flow  near  this  edge  in  more  detail,  introduce  local  (stretched)  coordinates 


x = ± r = 4, 

a4 ' a4 ' 


Z = X + iY  = — 
a4 


In  these  inner  variables  for  a —>  0 the  training  edge  of  the  plate  recedes  to  the  infinity  and  the  inner  flow 
becomes  as  shown  in  Fig.  6.  For  zero  cavitation  number  this  flow  is  completely  characterized  by  the  distance  s 
of  the  stagnation  point  from  the  leading  edge  of  the  plate.  The  inner  problem  is  solved  with  use  of  the  velocity 
hodograph  method  and  conformal  mappings. 

Introduce  the  hodograph  variable 

r^l_=ejEi0 

w1  q ' y ’ 

where  wt(Z,a ) is  the  conjugate  complex  velocity  (of  relative  fluid  motion)  in  the  inner  region  of  the  flow, 
q = \w‘\  and  9 = — arg'tF.  To  be  able  to  determine  the  flow  pattern  in  the  nearfield  using  hodograph  complex 
plane,  one  requires  additionally  an  estimation  of  the  behaviour  of  the  angle  9 of  the  tangent  to  the  (upper) 
boundary  of  the  cavity  at  the  downstream  infinity.  Essentially,  this  is  equivalent  to  estimating  of  the  one-term 
inner  limit  of  the  one-term  outer  expansion.  Replacing  the  outer  variables  in  the  expression  (61),  by  the  inner 
variables,  i.e.  x = a4X  and  y = a4Y,  and  expanding  for  X = 0(1)  and  a — > 0,  one  obtains 

y(x)  = y(a4X)  I =y°i(X,a)  = 

\ot—*0 


a4  [4^3/4  + |a2x5/4  _ _Lq4X7/4  + 


wherefrom  the  one-term  inner  expansion  of  the  cavity  contour  in  the  form 

Y = = [.V  ! + 0(a2)  (65) 

a4  3 

The  inclination  of  this  curve  for  X —>  oo  tends  to  zero,  so  that  one  can  assume  that  9 = 0 at  the  downstream 
infinity.  The  hodograph  plane,  shown  in  Fig.  7,  is  then  transformed  onto  an  auxiliary  plane  £ in  such  a way 
that  the  images  of  the  plate  and  the  cavity  be  found  on  the  real  axis  5I(  = £,  see  Fig.  7. 

. 1/  1\ 


To  complete  the  procedure  of  obtaining  the  inner  solution  it  is  necessary  to  relate  the  complex  potential  of 
relative  motion  of  the  fluid 

Fi  = ip  + iip,  (67) 

where  ip  and  'ip  are  correspondingly  the  velocity  potential  and  stream  function  (of  relative  motion)  with  the 
function  wl  by  means  of  an  intermediate  complex  plane  £.  The  complex  potential  plane  F*.  shown  in  Fig.  7,  is 
transformed  onto  (-plane  with  help  of  the  function 

Fi=(^V  (68) 

where  At  is  a real  constant.  Using  formulae  (63)  and  (66),  one  comes  to  the  following  relationship 

wi  = £ ± \J (2  - 1,  (69) 

where  the  plus  sign  is  used  for  |£|  < 1,  and  the  minus  sign  should  be  taken  for  j£j  > 1. 

The  relationship  between  Z and  w’  through  complex  variable  £ can  be  derived  in  the  following  way 


az  = z2-z1=  / --4. 


d Fi  d£ 
d£  wl 
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Using  formulae  (68)  and  (69)  in  combination  with  (70)  and  integrating  from  the  leading  edge  to  the  stagnation 
point,  one  finds 

12 

A = -a  (71) 

Making  use  of  (70)  one  can  find  the  following  relationship  between  Z and  ( = £ on  the  plate  and  the  cavity  in 
the  following  form 


• On  the  plate,  for  0 < X < a,Y=0-0 

X_24r  1 1 1(|2-1)3/2  3i 

7 ~ m2(£-  l)2  + i-  1 ~ 3 (f  - l)3  + 8J 

• On  the  plate,  for  a < X < oo.Y-0-0 

X 24  r 1 1 1(|2-1)3/2  ll 

a 17L2(|-1)2  £-1  3 (£-l)3  3J 

• On  the  cavity,  -1  < | < 1 

Z X Y 24  r 1 1 3 i (1  - £2)3/2 1 

a a + 1 a ~ 17  l.2(f  -l)2+£-l  + 8 + 3 (1  - £)3  J 


(72) 


(73) 


(74) 


The  latter  equation  provides  a parametric  relationship  between  Y = Y(£)  and  X = X(£)  on  the  cavity  contour 
in  inner  variables. 

Having  determined  the  relationship  Z = Z(Q,  described  by  (70),  with  account  of  (69)  and  (70)  one  can  find  an 
implicit  expression  for  distribution  of  the  velocity  along  the  plate  in  the  inner  region 


and 


X _ 24  r 2(,ui)2  2wl  l/l  + uV*  3i 

7 ~ 17  L (1  - -id)4  + (1  - 'id)2  + 3 V 1 ai)  +gj’ 

X__  24  r 2(ui)2  2 1 /I  +»1\3  li 

a 17  1.(1  — m*)4  (1  — u1)'2  3 V 1 — ul ) 3J  ’ 


for  0 < X < a,  Y = 0-0; 
for  a < X < oo,  Y = 0 — 0. 


(75) 


(76) 


3.3.  Matching  and  additive  composition 

To  determine  the  parameter  a,  entering  the  inner  solution,  one  uses  the  Van-Dyke  asymptotic  matching  principle, 
[2],  First  of  all,  find  a one-term  inner  representation  of  the  two-term  outer  solution2 

On  the  other  hand,  take  the  outer  representation  of  (75),  (76)  for  u1.  Introducing  x = a4X,  one  can  write 

~ 1 -cmi  + 0(t*2).  (78) 


Substituting  this  expression  into  (75)  and  (76)  expanding  it  for  a — > 0 and  x =fixed,  one  has 


1 — a 


/48a\!/4 

\17x) 


(79) 


Comparing  formulae  (77)  and  (79)  gives 


(80) 


wherefrom  one  can  see  that  the  distance  of  the  stagnation  point  from  the  leading  edge  is  equal  to  a4/ 4. 
Uniformly  valid  expression  for  the  velocity  uc  on  the  wetted  side  of  the  plate  can  be  found  by  means  of  the 
additive  composition  of  the  outer  (56)  and  inner  (75), (76)  solution,  The  latter  procedure  implies  adding  these 
solutions  and  substracting  their  common  part  (77),  i.e. 


+ u°  - u%°  = ul  + u°  - u°l  = ul  - 


a\fjx 
1 + — \fx 


(81) 


The  uniformly  valid  distribution  of  the  pressure  coefficient  is  calculated  by  means  of  the  formula 

C2P  = 1 - (wc)2  (82) 

Figures  8 and  9 illustrate  comparison  of  the  pressure  distributions,  calculated  using  the  formula  (82)  and  the 
exact  solution  for  a = 10°.  The  results  of  the  linear  theory  are  plotted  with  dot-and-dash  line.  The  upper  cavity 
shape  (in  stretched  coordinates)  in  the  immediate  vicinity  of  the  trailing  edge  of  the  supercavitating  plate  is 
shown  in  Fig.  10. 


2 Adding  preliminarily  the  velocity  of  the  incoming  flow 
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4.  Supercavitating  foil  with  a spoiler 


As  another  example  of  how  the  matched  asymptotics  technique  can  be  used  to  complement  the  outer  (linear) 
expansion  of  the  solution  in  a local  region  of  large  flow  perturbations  is  that  of  a flow  past  a supercavitating 
plate  with  a spoiler,  [5],  [6].  The  spoiler  is  represented  by  a plate  of  a small  relative  height  mounted  upon  the 
lifting  surface  either  at  the  trailing  edge  on  the  pressure  side  or  at  an  appropriate  station  on  the  suction  side. 
The  spoiler  is  usually  oriented  normally  to  the  oncoming  stream.  The  presence  of  the  spoiler  results  both  in  local 
pressure  rise  due  to  creation  of  the  stagnation  zone  in  its  vicinity,  and  redistribution  of  pressures  around  the 
entire  foil,  and,  eventually,  the  additional  lift  occurs.  Experiments  of  recent  years  have  shown  that  the  spoiler 
is  one  of  the  most  effective,  yet  relatively  simple  devicess  to  enhance  the  lifting  capacity  of  the  supercavitating 
hydrofoils  as  well  as  the  thrust  of  the  supercavitating  screw  propellers.  It  is  worth  mentionning  that  Professor 
Tulin  marked  spoilers  on  one  of  the  branches  of  his  famous  family  tree  of  supercavitating  flow  theory,  [6] . 

In  the  example  below,  considered  is  the  case  of  a two-dimensional  flow  past  a supercavitating  foil  with  a spoiler 
of  a small  relative  width  s at  a trailing  edge.  For  the  sake  of  illustrating  the  procedure,  the  simplest  case  of  zero 
cavitation  number  is  considered.  In  the  outer  region  (far  from  the  spoiler  in  terms  of  its  length)  appropriate 
linear  solutions  are  used  incorporating  an  admissible  (square  root)  singularity  of  unknown  strength  at  the 
trailing  edge.  In  the  inner  region  (in  the  vicinity  of  the  spoiler)  the  problem  is  reduced  to  that  of  a Kirelihoff 
type  separated  flow  past  a symmetric  wedge.  Asymptotic  matching  of  the  outer  and  inner  solutions  permits  to 
determine  hydrodynamic  characteristics  of  the  supercavitating  foil  with  the  spoiler. 

Let  the  supercavitating  foil  be  slightly  curved  and  oriented  to  the  flow  at  a small  angle  of  attack  a.  As  per 
foregoing,  assume  that  the  cavitation  number  a = 0.  The  spoiler  has  a small  relative  width  e « 1 and  is 
oriented  at  an  arbitrary  angle  ft  to  the  foil  at  the  trailing  edge. 

4-1-  Flow  near  the  spoiler  (inner  problem) 

In  the  region  near  the  spoiler  introduce  stretching  of  the  local  independent  variables 

X = — Y = - . Z = X + i Y.  (83) 

e e 

Assume  that  the  distance  from  the  flow  boundaries  is  of  the  order  of  the  chord,  i.e.  0(1).  Then,  the  pattern  of 
the  local  flow  does  not  depend  on  the  type  or  on  the  number  of  the  boundaries.  It  does  not  depend  either  on 
the  distance  of  the  trailing  edge  from  the  flow  boundaries,  and  represents  a flow  past  a semi-infinite  horisontal 
flat  plate  with  a spoiler  of  a unit  width.  Analytic  continuation  of  this  flow  into  the  upper  half-plane  leads  to  a 
problem  for  a separated  flow  of  unit  velocity  past  a symmetric  wedge.  The  solution  of  this  classical  problem  can 
be  obtained  by  the  methods  of  the  jet  theory  in  an  ideal  fluid.  Following  [7]  find  the  conjugate  complex  velocity 
wl{Z),  by  mapping  the  complex  potential  plane  Fl{Z ) = pl  + up13  and  the  logarithmic  velocity  hodograph 
plane  u>  = — hi  w = — In  \w\  + i 9 upon  auxiliary  plane  t (Fig.  11): 


F{Z)  = p0t2, 


2/3,  , 

UJ  = — ill 


(v1_^  + D’ 


(84) 


where  p0  is  a parameter  related  to  the  wedge  angle  ft  (in  radians)  by  means  of  the  following  relationship 

For  better  convergence  of  the  integral  in  the  denominator  it  is  practical  to  to  use  it  in  the  following  alternative 
form 

fjf rf+1)w%d.= 

/■  V-W.)  [(7n^+  if”  - sW'l  it  + 2W*"‘  (86) 

Jo  ' J ir-P/TT 

The  coordinates  of  the  free  boundary  of  the  cavity,  detaching  from  the  spoiler  (lower  cheek  of  the  wedge)  can 
be  calculated  with  use  of  the  formula 


Z = 


VT^¥  + l 

t 


\ 2/3/ vr 

J tdt  + exp(— i/3) 


(87) 


^of  relative  fluid  motion 
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In  order  to  perfom  the  matching  of  the  inner  solution  to  the  outer  solution  one  has  to  obtain  the  asymptotics 
of  the  inner  solution  far  the  spoiler.  With  the  purpose  to  perform  the  matching  of  the  ordinates  of  the  free 
boundary  of  the  cavity4  find  the  asymptotics  of  the  free  boundary  for  X —>  oo,  ( t — > — oo) 


where 


wherefrom 


Z(t)  — X(t)  + \Y (t)  ~ ‘2cp0  t2  + — t)  + C 

K/ 


C*  = 2 cp0 


1 i\2/5A  1 , 

J — 1 rdr 


Ti  T 


y 4 0 


Y = — ~ ( pQX , for  X — ► -boo 

£ 7T 


(88) 

(89) 

(90) 


Besides,  estimating  the  outer  limit  of  the  inner  solution,  one  can  show  that  the  outer  description  of  the  conjugate 
complex  velocity  has  a square-root  singularity  at  the  point  x = 1.  In  fact,  expanding  w(t)  and  Z(t)  for  large  t, 
one  can  obtain 

Z(t)~<p0t2  + 0(t), 


and 


(/  1 i\~  zp/ir  Vi  H 

Excluding  the  auxiliary  variable  t,  from  these  two  expressions,  one  can  find 


' dZ  ,yz 


On  the  foil 


Z = X = 


x - 1 


Finally,  one  can  evaluate  the  behaviour  of  w on  the  foil  near  the  trailing  edge 

W\ Xplyfe 


7T  y/ 1 - X 


for  X — > 1 — 0 


(91) 

(92) 

(93) 

(94) 

(95) 


4-2.  The  outer  ( linearized ) flow 

For  the  simplest  case  of  zero  cavitation  number  (a  = 0)  and  angle  of  attack  a « 1 the  flow  field  outside  of 
the  vicinity  of  the  spoiler  experiences  small  perturbations.  The  corresponding  linear  solution  for  the  conjugate 
complex  velocity  can  be  written  in  the  form 

»(C)  = 4-yS)__|i_  (96) 

Note  that  the  foregoing  expression  differs  from  the  solution  for  the  foil  without  spoiler  (42)  by  the  last  term, 
where  Bs  is  a real  constant.  This  latter  term  represents  an  admissible  solution  of  the  corresponding  Riemann- 
Hilbert  mixed  boundary  problem  which 

• has  zero  real  part  on  £ = £ > 1 and  C = £ < 0,  i.e.  on  the  boundary  of  the  cavity 

• has  zero  imaginary  part  on  1 > (,  = £ > 0,  i.e.  on  the  wetted  part  of  the  foil 

• yields  a square  root  singilarity  of  the  type  1 / flx  — 1 at  the  point  C = $ = 1—0,  corresponding  to  the 
point  x = 1 — 0.  In  other  words,  such  a solution  secures  matching  of  the  outer  solution  with  the  inner 
solution,  discussed  above 

In  fact  near  that  point  (£  — > 1 — 0,  and  x — > 1 — 0)  the  “spoiler”  term  behaves  as 
iBa  Bs  _ Bs  / 1 + yjx 

- 1)  x v .rv/l  v/.,-  y/v^v/T^ 

4 Note  that  the  matching  can  be  performed  with  respect  to  other  parameters  of  the  flow,  e.g 


v/2  Bs 
\/l  — a 


(97) 


pressure  or  velocity  on  the  foil 
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• yields  a (proper)  quarter-root  singularity  of  the  type  x x/4  on  the  wetted  side  of  the  foil  close  to  the 
leading  edge. 

In  fact,  near  the  latter  point  one  has 


iZL 


B, 


\M  - 1) 


'xsjl  - y/x 


Ba 

'xl/4 


(98) 


• complies  with  the  requirement  of  the  decay  of  the  perturbation  velocity  at  the  infinity 

The  real  constant  Bs  characterizes  the  strength  of  the  square-root  singularity  at  the  point  x = 1 and  has  to  be 
determined  from  the  matching  procedure.  Matching  of  the  linearized  outer  solution  and  the  inner  solution  can 
be  performed  in  terms  of  the  perturbation  velocities,  or  in  terms  of  the  ordinates  of  the  lower  free  boundary  of 
the  cavity  detaching  form  the  spoiler.  Using  the  former  option  and  employing  the  matching  principle,  one  should 
equate  the  outer  expansion  wl°  of  the  inner  solution  of  with  the  inner  expansion  wm  of  the  outer  solution5.  In 
other  words,  one  should  equate  the  expressions  (95)  and  (97).  As  a result,  the  constant  Bs  is  obtained  in  the 

form  

/?\/2 ^ 


B,  = 


~V~e 


(99) 


It  is  easy  to  verify  that  the  matching  could  have  been  done  in  terms  of  the  ordinates  of  the  lower  free  boundary 
of  the  cavity.  Find  the  asymptotics  of  the  outer  description  of  this  boundary6  near  the  spoiler,  i.e.  for  x — > 1 — 0. 
Using  the  kinematic  condition  on  the  cavity  u = 0 


Ay 


= v = — Shn  = 


Bs 


(£  = s/x.x  >!,£>!) 


d*  " ~ v^Cl  - 0 ' 

Integrating  this  expression  one  obtains 

fx  f ^ Ax  f v/7 

y = I ^ vAx  = j w(0^7  (1C  = 2 I v(|)£d£ 


Asymptotic  representation  of  (100)  for  x 
constant  B 


x-  1)  + ln[2v/x  — 1 + 2y  \fx{\/x  — 1)]}.  (100) 

1 — 0 is  with  use  of  the  previously  obtained  magnitude  of  the 


V 


— 2 V2BS  y/x^l  = 


(101) 


Re-writing  (90)  in  terms  of  the  outer  variable  x one  can  verify  complete  coincidence  of  the  resulting  expression 
with  (101).  This  fact  proves  correctness  of  the  matching  procedure. 

The  additional  lift  coefficient  due  to  presence  of  the  spoiler 


CVH  = -25? 


<ji  wa(z ) Az  = — J us(x)  dxr;  = 2Bit  = 2^/2 


The  drag  coefficient  due  to  presence  of  tli  spoiler  in  the  case  of  flat  plate  and  zero  angle  of  attack  is 

^oP2 


6’  = -CA 


<j)  w2(z)  Az  = — Q w2(C)^dC  = 27r B2 


(102) 


(103) 


Some  calculated  results  reflecting  dependence  of  the  spoiler  contributions  to  the  lift  and  drag  coefficients  on  the 
spoiler  installation  angle  fi  are  presented  in  Fig.  12.  For  the  case  of  non- zero  angle  of  attack  a ^ 0 the  resulting 
lift  and  drag  coefficients  can  be  determined  with  use  of  the  known  superposition  rules 


Cy  — CVa  + Cys 

C-x  = {\fCXa  + \/  CXs)2 

For  the  case  of  the  flat  plate  the  concrete  expressions  for  these  coefficients  become 

Cy  = — + 2 \J  2 ip„P\fe 


(104) 

(105) 

(106) 


5 In  the  overlap  regions 

6Based  on  the  possibility  of  linear  superposition,  it  is  sufficient  to  take  the  case  of  zero  incidence 
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Cx  = (\/|a  + 2./%Vi)2  (107) 

\ V l y 7 r 

The  spoiler  contribution  to  the  upper  boundary  of  the  cavity  can  be  calculated  by  integrating  the  corresponding 
“upwash”  = v for  z = x + iO,  Q = £,  £ = —y/x,  x > 0.  Using  (96),  one  has 


V = 


-vdF 

2 / t>(|)|d£ 

0 


2 B. 


[\Jx  + y/x  + ln( yjl  + yfx  - \[jx)\ 


(108) 


The  structure  of  asymptotic  solutions  can  be  used  for  some  useful  estimates.  One  can  evaluate  e.g.  optimal 
ratio  angle  of  attack  and  relative  length  of  the  spoiler.  Utilizing  expressions  (106)  and  (107)  one  can  write  the 
following  formula 


L r Cv  aaK  + bs 

dV£  = C~x  = (CccK  + de)2  + Cf/e 


(109) 


where  the  coefficients  aa,ca,be,de  can  be  easily  determined  comparing  (106)  and  (107).  Parameter  k is  defined 


as 


a 


(110) 


and  Cf  = Cf(Re)  is  a friction  coefficient,  which  can  be  calculated  as  a function  of  Reynolds  number  with  use  of 
an  appropriate  formula.  For  example,  assuming  fulll  turbulent  regime  of  the  flow  past  the  wetted  part  of  the 
cavitating  plate,  one  can  use  the  formula 


0.455 

Cf  ~ (log  Re)2  ™ 


(111) 


Differentiating  (109)  with  respect  to  k and  equating  the  result  to  zero,  one  obtains  the  following  optimal 
magnitude  of  k 


Kopt 


-Cgb€  + \J ( <‘a bE  - Onde)2  + Cf/e 
O' a (‘ce 


(112) 


Note  that  for  cj/e  — »■  0 one  can  assume  that 


Kopt 


d.  2vd>~-  Tj.i 

Ca  n 


(113) 


so  that  the  optimal  ratio  a/ \fe  is  negative.  The  latter  results  means  that  to  secure  maximum  lift-to-drag 
ratio,  one  has  to  provide  a negative  angle  of  attack,  for  example  for  a spoiler  normal  to  the  plate  (/ 3 = tt/2), 
atopt  = — 0.748\/e. 

The  value  of  Kopt  can  be  employed  to  calculate  the  maximum  lift-to-drag  ratio  of  a flat  plate  with  a spoiler  at  the 
trailing  edge.  However,  there  exits  a somewhat  contradicting  requirement  of  univalence  which,  in  simple  words, 
means  that  the  upper  side  of  the  cavity  generated  by  the  optimal  combination  of  the  spoiler  and  (negative)  angle 
of  attack  should  not  intersect  the  plate7.  One  can  find  the  domain  of  univalence  by  considering  the  thickness  of 
the  cavity  with  respect  to  the  plate.  Using,  the  corresponding  expressions  (61)  and  (108)  for  the  contributions 
to  the  ordinates  of  the  upper  side  of  the  cavity  due  to  angle  of  attack  and  presence  of  spoiler,  one  finds  the 
aforementionned  cavity  thickness  in  the  form 


Vt(x)  - ya  + ys  - afa(x)  + ^ y/efe(x) 

7 r 

(114) 

where 

fa(x)  = -[(1  + 2\fx)\J x + \fx  - ln(^l  + Vx+  y/v^)] 

(115) 

fe{ x)  = 2 [\Jx  + y/x  + ln( -y/ 1 + y/x  ~ \[jc)} 

(116) 

Requiring  that  the  cavity  thickness  be  non-negative  at  the  trailing  edge  of  the  plate  (yt(  1)  < 0),  one  obtains 
the  expression  for  the  domain  of  univalence  in  terms  of  of  the  relationship  between  angle  of  attack  and  length 
of  the  spoiler 

K = 4=  < = -NoPyfcffi,  (117) 

/«(!) 

' Condition  of  univalence  for  this  case  was  first  found  by  A.S.  Achkinadze  and  G.M.  Fridman 
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where 


4\/2  \/2  + ln(\/2-l) 
7T  3\/2  - ln(\/2  + 1) 


0.285 


(118) 


The  optimal  and  “univalent” ratios  a/s/e  are  plotted  versus  the  spoiler  installation  angle  0 in  Fig.  13.  This 
Figure  clearly  shows  that  in  a certain  range  of  0 the  optimal  ratio  of  the  angle  of  attack  and  the  length  of 
the  spoiler  cannot  be  achieved  in  reality.  The  ( “pedestrian” ) maximum  of  lift-to-drag  ratio  derived  through 
optimisation  of  the  ratio  a/  yje  is  plotted  in  Fig.  14  together  with  the  lift-to-drag  ratio  attainable  within  the 
regime  of  univalence. 
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Figures  to  “Supercavitating  flows:  small 
perturbation  theory  and  matched  asymptotics” 


Fig.  1 Selected  nonlinear  cavity  closure  models  (1  - Riaboushinsky  model.  2 - Efros-Gilbarg  model,  3 - Tulin 

double-spiral  vortex  model,  4 - Wu-Fabula  model) 


Fig.  2 Linearized  “physical”  complex  plane  z and  auxiliary  complex  plane  ( for  the  case  of  closed  cavity  with 

singular  termination,  the  ground 


Lift  Coefficient  for  Different  Closure  Schemes 
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Fig.  3 Linearized  “physical”  complex  plane  z for  the  case  of  smooth  cavity-wake  transition  (Wu-Fabula  cavity 

closure  scheme) 


Fig.  4 Linearized  “physical”  complex  plane  z and  auxiliary  complex  plane  £ for  the  case  of  cavity  termination 

on  Tulin  double-spiral  vortices 
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Fig.  5 Lift  coefficient  ( Cy/a ) of  a supercavitating  flat  plate  versus  relative  cavitation  number  cr/a  for 
different  linear  analogues  of  the  cavity  closure  models  (1  - Tulin  double-spiral  vortex  scheme,  2 - closed  cavity 

scheme,  3 - Wu-Fabula  scheme) 


Pressure  Coefficient  ~ w Pressure  Coefficient 
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Abcissa 


ig.  8 Pressure  distribution  along  supercavitating  fiat  plate  (a°  = 10°)  (solid  line  - uniformly  valid  MAE, 
dashed  line  - exact  solution,  dot-and-dash  line  - linear  theory) 


Abcissa  Measured  from  Leading  Edge 


Fig.  9 Pressure  distribution  in  the  vicinity  of  the  leading  edge  of  a supercavitating  flat  plate  for  a°  = 10' 

(solid  line  - MAE,  dashed  line  - exact  solution) 


Fig.  10  Upper  cavity  shape  in  stretched  coordinates  in  the  immediate  vicinity  of  the  leading  edge  of 

supercavitating  flat,  plate 
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Fig.  11  “Physical”  plane,  complex  potential  plane  and  hodograph  plane  for  the  flow  problem  solution  near 

the  spoiler 
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Fig.  12  Lift  and  drag  coefficients  of  a supercavitating  fiat  plate  with  a spoiler  versus  spoiler  installation  angle 
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Fig.  13  Optimal  and  univalent  ratios  a/ \fe  versus  the  spoiler  installation  angle  for  a supercavitating  flat 

plate  with  a spoiler  at  the  trailing  edge 


Maximum  Lift-to-Drag  Ratio 
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Fig.  14  Maximum  “pedestrian”  and  “univalent”  lift-to-drag  ratios  of  a supercavitating  flat  plate  with  a 

spoiler  versus  the  spoiler  installation  angle 


